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Abstract 


Star-forming processes strongly influence the ISM chemistry. Nowadays, many high-quality databases are 
available at millimeter wavelengths. Using them, it is possible to carry out studies that review and deepen previous 
results. If these studies involve large samples of sources, it is preferred to use direct tools to study the molecular 
gas. With the aim of testing these tools such as the use of the HCN/HNC ratio as a thermometer, and the use of 
НСО”, HC4N, N;H* and СН as “chemical clocks,” we present a molecular line study toward 55 sources 
representing massive young stellar objects at different evolutionary stages: infrared dark clouds (IRDCs), high- 
mass protostellar objects (HMPOs), hot molecular cores (HMCs) and ultracompact H II regions. We found that the 
use of the HCN/HNC ratio as a universal thermometer in the ISM should be taken with care because the HCN 
optical depth is a big issue that can affect the method. Hence, this tool should be utilized only after a careful 
analysis of the HCN spectrum, checking that no line, neither the main nor the hyperfine ones, presents absorption 
features. We point out that the analysis of the emission of H'CO*, HC3N, NH * and СН could be useful to trace 
and distinguish regions among IRDCs, HMPOs and HMCs. The molecular line widths of these four species 
increase from the IRDC to the HMC stage, which can be a consequence of the gas dynamics related to the star- 
forming processes taking place in the molecular clumps. Our results not only contribute with more statistics, acting 
as a probe of such chemical tools, useful to obtain information in large samples of sources, but also complement 
previous works through the analysis of other types of sources. 
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1. Introduction 


The sites in which the stars form are characterized by a rich 
and complex chemistry. The smallest gaseous fragments within 
a molecular cloud, known as hot molecular cores (HMCs), 
which are related to the formation of massive stars, are the 
chemically richest regions in the interstellar medium (ISM) 
(e.g., Herbst & van Dishoeck 2009; Bonfand et al. 2019). 
Molecules and chemistry are ubiquitous along all the stages 
that a forming star goes through, and moreover, the star- 
forming processes strongly influence the chemistry of such 
environments (Jørgensen et al. 2020). For instance, as material 
collapses and becomes ionized by the young massive stars and 
shocked by jets and outflows, temperature and density can 
change drastically, leading to the formation and destruction of 
molecular species. Thus, observing molecular lines and study- 
ing their emission and chemistry are important to shed light on 
the different stages of star formation and to characterize 
physical and chemical conditions. 

To explore the star-forming processes in different environ- 
ments, one should work with accurate values of the physical 
properties, such as the case of the kinetic temperature (Тк). For 
example, calculating the excitation temperature (Tex) of carbon 
monoxide '^C'60, the Ty of a molecular cloud can be roughly 


estimated if it is assumed that there is a complete thermalization 
of the lines (То, = Тк). Ammonia (NH3) has been found in 
different interstellar environments: from dark quiescent clouds, 
circumstellar envelopes, and early stages of high and low 
]uminosity star formation to planetary atmospheres and external 
galaxies (Betz et al. 1979; Martin & Ho 1979; Ho & 
Barrett 1980; Rathborne et al. 2008; Takano et al. 2013; 
Bógner et al. 2022). Estimating the NH; rotational temperature 
usually results in a reliable indicator of the kinetic temperature. 

Hydrogen cyanide (HNC) and isocyanide (HNC) are two of 
the most simple molecules in the ISM, first detected almost 
fifty years ago (Snyder & Buhl 1971, 1972). These isomers 
have a linked chemistry, and differences in the spatial 
distributions in which they lie within a cloud can reflect the 
gas chemical conditions and the evolution of the star-forming 
regions (Schilke et al. 1992). Recently, Hacar et al. (2020) 
proposed the HCN-to-HNC integrated intensity ratio as a direct 
and efficient thermometer of the ISM with an optimal working 
range 15 K S Тк < 40 K. The authors performed an analysis of 
such isomers throughout the Integral Shape Filament in Orion, 
deriving an empirical correlation between the HCN/HNC ratio 
and the kinetic temperature Тк. Based on the analysis of such 
correlation toward many dense molecular clumps from the 
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MALTO0 survey, they proposed that the HCN/HNC thermo- 
meter can be extrapolated for the analysis of the ISM in 
general, particularly in star-forming sites, aiming to explore it 
toward different regions and sources along the ISM. 

Molecular species that emerge and are destroyed during the 
birth of stars can be used to track the star-forming processes 
within molecular clumps and cores (Stephens et al. 2015; 
Urquhart et al. 2019). Comparisons between column densities 
and molecular abundance ratios that can be used to estimate the 
age and mark the evolutionary stages of star-forming regions are 
known as "chemical clocks." As Sanhueza et al. (2012) pointed 
out, only molecules that exhibit differential abundances with 
time can be used to evaluate the evolutionary status of a star- 
forming region. In general, chemical clocks have been studied in 
depth in low-mass star-forming regions, but this effort has been 
less developed in the context of high-mass star-forming regions. 
As shown by Yu & Wang (2015), in the context of an analysis of 
chemical clocks it is important to not only study abundance 
ratios of such molecules, but also the integrated line intensities, 
the line widths, among other parameters. For instance, the 
molecular line widths (Av WEM) are related to the gas 
kinematics of the molecular clump interior, regarding turbulence, 
outflows and shocks among other processes (Sanhueza et al. 
2012) which can give information about the evolutionary stage 
of a massive young stellar object (MYSO). 

Some interesting molecular species for probing physical and 
chemical properties of star-forming regions are diazenylium 
(N>H*) and ethynyl radical (C,H). Both molecules seem to be 
good tracers of dense gas in the early stages of the star-forming 
evolution (Beuther et al. 2008; Sanhueza et al. 2012, 2013), 
N>H* traces cold gas due to its resistance to depletion at low 
temperatures (Li et al. 2019), and the latter, additionally, can 
indicate the presence of a photodissociation region (PDR), 
where ultraviolet (UV) photons from young and hot massive 
stars irradiate acetylene to produce C5H (Fuente et al. 1993; 
Nagy et al. 2015; García-Burillo et al. 2017). НСО“ and 
Hcot (formylium species) are usually employed to 
investigate infall motions and outflow activity (Rawlings 
et al. 2004; Veena et al. 2018), and HC3N is helpful to explore 
gas associated with HMCs (Bergin et al. 1996; Taniguchi et al. 
2016; Duronea et al. 2019). 

The mentioned molecules are among the brightest lines, and 
they were called molecular fingerprints in a study of a large 
sample of molecular clumps (Urquhart et al. 2019). As the 
authors pointed out, such molecules are able to trace a large 
range of physical conditions including cold and dense gas 
(HNC, H"CO*, HCN, HN'*C, НСМ), outflows (HCO*), 
early chemistry (СН), and gas associated with protostars and 
young stellar objects (YSOs) (HC3N, and cyclic molecules). 
Thus, the analysis of such molecules gives us a significant 
amount of scope to search for differences in the chemistry as a 
function of the evolutionary stage of the star formation taking 
place within molecular clumps. For instance, Yu & Wang 
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(2015) studied 31 extended green object (EGO) clumps with 
data from the MALT90 aiming to better understand the 
chemical processes that take place in the evolution of massive 
star formation. They classified the sources and made a 
molecular line study over 20 MYSOs and 11 HN regions. 
Through the comparison of integrated intensities, line widths 
and column densities, derived from the emission of МН? and 
C,H with those of the formylium (НСО?) апа cyanoacety- 
lene (HC3N), they suggested that NH” and CoH could act as 
efficient chemical clocks. They found that the N3H* and СН 
column densities decrease from MYSOS to H П regions, and the 
[N;H*]/[H^CO^] and [C;H]/[HPCO *] abundance ratios 
also decrease with the evolutionary stage of the EGO clumps. 
In addition they ascertained that the velocity widths of МН”, 
C,H, НЕСО" and HC4N are comparable to each other in 
MYSOS, while in Н II regions the velocity widths of N;H* and 
СН tend to be narrower than those of НСО? and HC3N. 

Nowadays, many high-quality databases are available that 
were generated from observations obtained with the most 
important (sub)millimeter telescopes such as the IRAM 30m 
Telescope and the Atacama Large Millimeter/submillimeter 
Array (ALMA), among others. For instance, using this kind of 
data, it is possible to carry out new chemical studies that, in 
turn, review and deepen previous results. If these studies 
involve large samples of sources, it is preferred to use direct 
tools to study and probe the molecular gas like what is 
presented in Hacar et al. (2020) and Yu & Wang (2015). 

Infrared dark clouds (IRDCs) are massive, dense and cold 
clumps that may harbor budding stars, while high mass 
protostellar objects (HMPOs) are already protostars accreting 
material from their surroundings; at this stage, both temperature 
and density increase, but it is thought that at the beginning, they 
are chemically poor from an evolutionary point of view. HMCs 
are considered hotter sources where the chemistry is prolific as a 
consequence of embedded and evolved HMPOs (Giannetti et al. 
2017) until they eventually reach the last stage here considered: 
ultracompact H п (UCHN) regions. The stars responsible for the 
UCHII regions have generally finished their accretion process 
and have begun to ionize the gas around them. All the processes 
involved in this evolutionary path impact the molecular gas that 
eventually can be investigated through the emission of molecular 
lines. It is important to highlight that the mentioned phases in the 
massive star formation may not have well-defined limits, and 
sometimes one determined source may have some physical 
conditions overlapping with those of another kind of source 
(e.g., Beuther et al. 2007; Gerner et al. 2014). Moreover, 
different kinds of sources may be embedded within the same 
molecular clump. Given this complex scenario, comparative 
studies with many sources are needed to analyze the involved 
physics and chemistry. 

With the aim of testing the presented tools concerning the 
use of the HCN/HNC ratio as a thermometer (following the 
methodology published in Hacar et al. 2020), and the analysis 
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Table 1 
Analyzed Molecular Lines 


Molecule Transition Rest Frequency (GHz) 
HP"co* 1-0 86.7542 
C-H ° 1-0 3/2-1/2 F = 2-1 87.3169 
HCN 1-0 F = 1-1 88.6304 
1-0 F = 2-1 88.6318 
1-0 F = 0-1 88.6339 
HNC 1-0 90.6635 
HC3N 10-9 90.9790 
NH 1-0 F; = 2-1 F = 3-2 93.1737 
Note. 


“ Species with several observed hyperfine lines. Only the used line of each 
molecule is included in this work. 


of НСО“, HC3N, N;H* and СН (following the methodol- 
ogy reported in Yu & Wang 2015), we present this study 
toward 55 sources representing MYSOs at different evolu- 
tionary stages as described above. After the presentation of the 
data and the analyzed sources (Section 2), the paper is 
structured as follows: a reporting of the results regarding each 
tested tool (Section 3), their respective discussion (Section 4) 
and a summary of the main results (Section 5). 


2. Data and Analyzed Sources 


We made use of the catalog “IRAM 30 m reduced spectra of 
59 sources" J/A+A/563/A97 in the VizieR database! (Gerner 
et al. 2014), from which the molecular data were obtained and 
the sources were selected to perform our study. The particular 
molecular lines analyzed in our work are presented in Table 1. 
These data, obtained by Gerner et al. (2014), with the 30m 
IRAM telescope, are single pointing spectra observed toward 
each source. The data used here correspond to the 86-94 GHz 
band, which have an angular and spectral resolution of 29" and 
0.6 km 5! respectively, and typical 1с rms values of about 
0.03 K. 

The sources, classified as IRDCs, HMPOs, HMCs and UCHII, 
are presented in Table 2, where the coordinates and distances are 
included. The classification is based on the guidelines introduced 
by Beuther et al. (2007) and Zinnecker & Yorke (2007) 
according to the physical conditions of the evolutionary sequence 
in the formation of high-mass stars. IRDCs are sources that 
consist of cold and dense gas and dust that emit mainly at (sub) 
millimeter wavelengths. According to Gerner et al. (2014), the 
IRDCS of this sample consist of starless IRDCs as well as IRDCs 
already starting to harbor point sources at jum-wavelengths. 
HMPOs are sources hosting an actively accreting massive 
protostar(s), which shows an internal emission source at mid- 
infrared. HMCs are much warmer than HMPOs and can be 
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distinguished from a chemical point of view. At this stage, the 
central source(s) heats the surroundings evaporating molecular- 
rich ices. Finally, the UV radiation from the embedded massive 
protostar(s) ionizes the surrounding gas and gives rise to an 
UCHI region. Even though there might be overlaps among 
HMPO, HMC and even the UCHII stage, such classification was 
followed by Gerner et al. (2014) based on physical quantities, 
especially the temperature, which rises from IRDCs to HMPOs to 
UCHII regions, and based on the chemistry that can differentiate 
HMPOs (chemical poorer sources) from HMCs (molecularly 
richer chemistry sources). The sample includes 19 IRDCs, 20 
HMPOs, 7 HMCs and 9 UCHII regions. In addition to the 
possible overlaps among the types of sources, it is possible that 
within the beam of the observational data lies more than one type 
of source. This is the case for HMC034.26, a region hosting 
several HII regions and a hot core. Thus, it is necessary to be 
careful with any analysis in relation to the evolutionary stages. 

The quality of these IRAM data and the variety of sources in 
the sample included in such a catalog allow us to perform a 
new analysis of several molecular species in order to probe the 
molecular gas conditions and the chemistry related to IRDCs, 
HMPOs, HMCs and ЏСНП regions with the aim of probing 
chemical tools. 

Additionally, we used infrared (IR) data to complement the 
information obtained from the molecular lines with the aim of 
probing, in these wavelengths, the activity of the star-forming 
regions. Data at the Ks band extracted from the UKIRT 
Infrared Deep Sky Survey (UKIDSS) (Lawrence et al. 2007) 
and IRAC-Spitzer 4.5 jum data obtained from the Galactic 
Legacy Infrared Midplane Survey Extraordinaire (GLIMPSE) 
survey (Churchwell et al. 2009) were used. 


3. Results 
3.1. Infrared and Submillimeter Continuum Emission 


Given that jets and outflows affect the star-forming regions’ 
chemistry, we looked for evidence of such processes by using 
IR data to have additional information on each source, besides 
its classification, to interpret and complement the molecular 
analysis. Data at the Ks band extracted from the UKIDSS were 
used to search for signs of jets (e.g., Paron et al. 2022) mainly 
in HMPOs, and also in HMCs and ЏСНП regions. IRAC- 
Spitzer 4.5 jum data, obtained from the GLIMPSE survey, were 
employed to analyze whether the sources present extended 
4.5 um, likely indicating outflow activity (e.g., Davis et al. 
2007). Additionally, each source was checked if it is cataloged 
as an “EGO—likely MYSO outflow candidate" according to 
the catalog of Cyganowski et al. (2008). 

In Table 2 (Cols. 6 and 7), we indicate the kind of emission at 
the Ks and 4.5 uum bands observed within the beam size of the 
IRAM data centered at the position of each source. If one or 
several point sources appear within the IRAM beam, it is 
indicated as "points"; extended emission, suggesting the presence 
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Table 2 


Sample of Analyzed Sources 
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Name R.A. (J2000) Decl. (J2000) Distance (kpc) Type Ks 4.5 um Fpeak(870 um) (Ју Беат 1) 
IRDCOI 1.1 18 10 28.4 —19 22 34 3.60 IRDC points point 1.76 
IRDCO028.1 18 42 50.3 —04 03 20 4.80 IRDC points no E 
IRDC028.2 18 42 52.1 —03 59 54 4.80 IRDC points ext. 

IRDC048.6 19 21 44.4 13 49 24 2.50 IRDC points no 

IRDCO079.1 20 32 22.0 40 20 10 1.00 IRDC points points 

IRDC079.3 20 31 57.7 40 18 26 1.00 IRDC no points 

IRDC18151 18 17 50.3 --12 07 54 3.00 IRDC ext. no E 
IRDC18182 18 21 15.0 —14 33 03 3.60 IRDC points points 0.99 
IRDC18223 18 25 08.3 --12 45 27 3.70 IRDC points no 1.58 
IRDC18306 18 33 32.1 —08 32 28 3.80 IRDC points no 0.68 
IRDC18308 18 33 34.3 —08 38 42 4.90 IRDC points no E 
IRDC18310 18 33 39.5 —08 21 10 5.20 IRDC points no 0.92 
IRDC18337 18 36 18.2 —07 41 00 4.00 IRDC points no 0.66 
IRDC18385 18 41 17.4 —05 09 56 3.30 IRDC points no 0.53 
IRDC18437 18 46 21.8 —02 12 21 7.30 IRDC points points 0.68 
IRDC184.1 18 48 02.1 —01 53 56 6.40 IRDC points no 1.06 
IRDC184.3 18 47 55.8 —01 53 34 6.00 IRDC points no 2.06 
IRDC19175 19 19 50.7 14 01 23 1.10 IRDC points no 0.48 
IRDC2008 1 20 10 13.0 27 28 18 0.70 IRDC no x E 
HMPO18089 18 11 51.6 —17 31 29 3.60 HMPO ext. points+ext. 9.25 
HMPO18102 18 13 11.3 —18 00 03 2.70 HMPO points points 3.66 
HMPO18151 18 17 58.1 —12 07 26 3.00 HMPO jets points+ext. EE 
HMPO18182 18 21 09.2 --14 31 50 4.50 НМРО diff. points+ext. - EGO 5.66 
HMPO18247 18 27 31.7 —11 45 56 6.70 HMPO points points+ext. 1.87 
HMPO18264 18 29 14.6 —11 50 22 3.50 HMPO ext. ext. - EGO 7.43 
HMPO18310 18 33 48.1 —08 23 50 5.20 HMPO points points 2.12 
HMPO18488 18 51 25.6 00 04 07 5.40 HMPO points points 2.80 
HMPO18517 18 54 14.4 04 41 40 2.90 HMPO points ext. 5.26 
HMPO18566 18 59 10.1 04 12 14 6.70 HMPO diff. points+ext. 3.83 
HMPO19217 19 23 58.8 16 57 44 10.50 HMPO diff. points+ext. 3.06 
HMPO19410 19 43 11.0 23 44 10 2.10 HMPO jets ext. 5.31 
HMPO20126 20 14 26.0 41 13 32 1.70 HMPO jets x 

HMPO20216 20 23 23.8 41 17 40 1.70 HMPO jets ext. 

HMPO20293 20 31 12.9 40 03 20 1.30 HMPO diff. ext. 

HMPO22134 22 15 09.1 58 49 09 2.60 HMPO diff. ext. 

НМРО23033 23 05 25.7 60 08 08 3.50 НМРО x x 

HMPO23139 23 16 10.5 59 55 28 4.80 HMPO x x 

HMPO23151 23 17 21.0 59 28 49 5.70 HMPO x x 

НМРО23545 23 57 06.1 65 24 48 0.80 НМРО x x E 
HMC009.62 18 06 15.2 —20 31 37 5.70 HMC diff. ext. 12.48 
HMC010.47 18 08 38.2 --19 51 50 5.80 HMC points points 35.01 
HMC029.96 18 46 04.0 —02 39 21 7.40 HMC points ext. 12.01 
HMC031.41 18 47 34.2 —01 12 45 7.90 HMC points ext. 61.68 
HMC034.26 18 53 18.5 01 14 58 4.00 HMC diff. ext. 55.65 
HMC045.47 19 14 25.7 11 09 26 6.00 HMC ext. ext. - EGO 5.49 
НМСО75.78 20 21 44.1 37 26 40 4.10 НМС diff.+ext. ext. e 
UCH005.89 18 00 30.4 —24 04 00 2.50 UCHII ext. ext. 25.73 
UCH010.10 18 05 13.1 —19 50 35 4.40 UCHII ext. points 0.92 
UCH010.30 18 08 55.8 —20 05 55 6.00 UCHI diff. ext. - EGO 7.67 
UCH012.21 18 12 39.7 —18 24 20 13.50 UCHI points points 11.58 
UCH013.87 18 14 35.8 —16 45 43 4.40 UCHN diff.+ext. ext. 5.52 
UCH030.54 18 46 59.3 —02 07 24 6.10 UCHI points points+ext. 2.24 
UCH035.20 19 01 46.4 01 13 25 3.20 UCHI points+diff. x e 
UCH045.12 19 13 27.8 10 53 37 6.90 UCHI jets ext. 7.58 
UCH045.45 19 14 21.3 11 09 14 6.00 UCHI ext. ext. 3.90 
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Figure 1. Near-IR Ks emission obtained from the UKIDSS database (top panels) and IRAC-Spitzer 4.5 шт emission (bottom panels) toward HMPO 20126 and HMC 
45.47, on the left and right respectively. In all cases, the dashed circle represents the position and the beam size (29") of the IRAM observations. 


of jets and outflows, is signified as "ext." and if this extended 
emission presents a jet-like morphology, it is stated as “jets.” 
Diffuse emission without any clear morphology is indicated as 
"diff" While “no” means that there is no emission, an “x” 
indicates there are no data at the corresponding wavelength. In 
the case of the 4.5 um emission, if the source is cataloged as an 
“EGO—likely MYSO outflow candidate,” this is indicated with 
“EGO” in Col. 7. As an example, Figure | displays two sources 
in the Ks band: on the left HMC45.47 shows a jet-like feature 
and some extended emission at the 4.5 шт band, while on the 
right HMC 45.47 appears with extended emission at both bands, 
and it is cataloged as an EGO. 


Additionally, in Table 2 (Col. 8), we include the peak flux of 
the submillimeter emission obtained from the APEX Telescope 
Large Area Survey of the Galaxy (ATLASGAL) compact 
source catalog (Contreras et al. 2013; Urquhart et al. 2014) for 
sources that lie within the IRAM beam centered at the 
coordinates indicated in Cols. 2 and 3. 


3.2. HCN/HNC Ratio and Tx 


Following the same procedure as done by Hacar et al. (2020) 
when they test the HCN/HNC kinetic temperature in a large 
sample of dense molecular clumps extracted from the MALT90 
survey (see their Section 4.4), we derived the Tx for each source 
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of the sample presented here. The correlation to derive Tk 
proposed by the authors is as follows: if the integrated intensity 
ratio of the isomers HCN and HNC (hereafter ICN) /I(HNC)) 
is <4, Ty is estimated as Тк = 10 x I(HCN)/I(HNC); and if 
I(HCN)/I(HNC) > 4, Tk =3 x (K(HCN)/I(HNC) —4) + 40 is 
used. As done by Hacar et al. (2020), we integrated the HCN 
(J = 1-0), including all hyperfine components and HNC 
(J = 1-0) over a given velocity interval (in Appendix we include 
some spectra indicating the integrated area in both isomers). The 
obtained integrated emission values, and the respective ratios, 
are presented in Table 3 (Cols. 2, 3 and 4, respectively), and the 
obtained kinetic temperatures following the mentioned relations 
are presented in Table 4. Typical errors in I(HCN) and I(HNC) 
are about 0.2 and 0.1 K km 5 ", respectively, which yield errors 
between 1 and 5%, and in some cases at most 10%, in the Tx. 
From now on, the obtained errors in the measured and derived 
parameters appear as bars in the figures, and for a better display, 
we do not include them in the tables. 

In the further analysis and comparisons, following Hacar 
et al. (2020), we do not consider values of kinetic temperatures 
obtained from the HCN/HNC ratio (hereafter T,(HCN/HNC)) 
lower than 15 K. These values are due to I(HCN)/I(HNC) 
ratios close to the unity, which according to the authors, the 
uncertainties in the method grow for such cases, likely due to 
the combination of excitation and opacity effects. In 
Section 4.1.1, based on the inspection of the HCN and HNC 
spectra of sources for which T,(HCN/HNC) < 15 К, we 
discuss this issue. 

In order to compare the Tk (HCN/HNC), we sought the dust 
temperature of each source. Dust temperature (Та) values 
were obtained from the maps? generated by the point process 
mapping (PPMAP) algorithm (Marsh et al. 2015) done to the 
Herschel infrared Galactic Plane Survey (Hi-GAL) maps in the 
wavelength range 70—500 шт (Marsh et al. 2017). Taust values 
are included in Table 4, and they represent average values on 
the dust temperature map over the IRAM beam size centered at 
the source position. Additionally, we also compared with the 
ammonia kinetic temperature (Tk(NH3)) derived by Urquhart 
et al. (2011) in several sources of the sample analyzed here that 
are contained in the Red MSX Source survey (Hoare et al. 
2005; Mottram et al. 2006; Urquhart et al. 2008). 

Figure 2 exhibits the comparisons between Tk(NH3) and 
Taust (upper panel), Taust and Tk(HCN/HNC) (middle panel), 
and Tx(NH3) and Tx(HCN/HNC) (bottom panel). Table 5 
presents the average temperature values with their errors. As 
mentioned above, sources with values T«(HCN/HNC) < 15 К 
were not considered, either in the figures or in the calculated 
average. 


2 http: //www.astro.cardiff.ac.uk /research/ ViaLactea / 
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3.3. Molecules as “Chemical Clocks"? 


Following Yu & Wang (2015), we performed a similar 
analysis using the HNC, CH, HC3N, НСО? and N>H* 
emissions measured toward the sources presented in Table 2. 
This is: (1) we compared the molecular emissions with the flux 
at the submillimeter continuum, (2) we analyzed the line widths 
of the C,H, HC3N, H'CO* and N-H*, and (3) we analyzed 
the column densities (extracted from Gerner et al. 2014) of CH 
and NoH* by comparing with that of the НСО“. 

Figure 3 displays the integrated intensities of HNC, СН, 
HC3N and N>H* (presented in Table 3) versus the 870 шт 
peak flux obtained from the ATLASGAL compact source 
catalog (presented in Col. 8 in Table 2). 

We fitted the line emission of СН, NoH*, НСО? and 
HC3N with Gaussian functions to obtain the full width at half 
maximum (FWHM) line widths (Av) in each source (in 
Appendix we include some spectra showing the Gaussian 
fittings). In the case of the СН and N;H^, which have 
hyperfine components, the fitting was done with multiple 
Gaussian functions, and the Av values used in this analysis 
correspond to the main components. These values are presented 
in Table 6. The absence of values in some sources is due to a 
lack of emission or that the hyperfine components are 
completely blended with the main component. Table 7 
expresses the average line widths for each kind of source, 
and with the aim of analyzing their behavior, following the 
analysis reported in Yu & Wang (2015), we present plots that 
display relations among the measured Av in Figure 4. 

From Gerner et al. (2014) we obtained the column densities 
of CoH, МН” and НСО? for each source. We used the 
column densities obtained as "iteration 1" according to the 
model used by the authors, which are the values derived with 
the mean temperatures from the best-fit models of "iteration 0." 
Values indicated as upper limits are not included in our 
analysis. Figure 5 displays the column densities of C,H and 
N5H* versus the H^CO" column density. Error bars are not 
included, given that the authors did not inform them in their 
work. Finally, Figure 6 depicts the relative abundance of 
[C>H]/[H'*CO*] versus [N)2H *]/[HP?CO *]. 


4. Discussion 


Unlike low-mass stars, the formation of high-mass stars is 
not yet fully understood, and different scenarios have been 
proposed (see e.g., Tan et al. 2014; Motte et al. 2018). Probing 
the chemical conditions of high-mass star-forming regions at 
different evolutionary stages is an important issue to advance in 
knowledge on the formation of this kind of star, and to explore 
such chemical conditions, it is necessary to use efficient tools. 

In this context, we decided to test the HCN/HNC ratio as a 
thermometer of the gas, as well as the use of the emission of 
H'CO*, HC4N, МН” and CH as "chemical clocks.” The 
obtained results from each study are discussed in what follows. 
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Integrated line Intensities (Units in K km s7!) and the Ratio between HCN and HNC 


Table 3 
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Source 


IRDCO11.1 
IRDCO28.1 
IRDCO28.2 
IRDC048.6 
IRDCO79. 1 
IRDCO079.3 
IRDC18151 
IRDC18182 
IRDC18223 
IRDC18306 
IRDC18308 
IRDC18310 
IRDC18337 
IRDC18385 
IRDC18437 
IRDC184.1 
IRDC184.3 
IRDC19175 
IRDC20081 
HMPO18089 
HMPO18102 
HMPO18151 
HMPO18182 
HMPO18247 
HMPO18264 
HMPO18310 
HMPO18488 
HMPO18517 
HMPO18566 
HMPO19217 
HMPO19410 
HMPO20126 
НМРО20216 
НМРО20293 
НМРО22134 
НМРО23033 
НМРО23139 
HMPO23151 
HMPO23545 
HMC009.62 
HMC010.47 
HMC029.96 
HMC031.41 
HMC034.26 
HMC045.47 
HMC075.78 
CH005.89 
CH010.10 
CH010.30 
CH012.21 
CH013.87 
CH030.54 
CH035.20 
CH045.12 
CH045.45 


сесесесеесоа 


Notes. 


KHCN) 


6.83 
16.70 
8.71 
3.36 
7.96 
6.40 
49.10 
6.03 
18.90 
1.28 
14.30 
2.59 
1.08 
6.51 
6.93 
0.60 
11.50 
2.46 
6.62 
10.40 
22.60 
49.10 
34.00 
10.60 
74.10 
8.60 
7.31 
48.20 
11.20 
18.30 
38.70 
55.70 
15.90 
29.70 
13.30 
35.70 
31.10 
14.50 
6.79 
45.10 
110.30 
49.20 
19.70 
24.70 
23.20 
33.90 
130.50 
4.48 
95.10 
13.60 
31.30 
21.20 
26.40 
43.40 
22.10 


KHNO) 


7.27 
12.50 
8.34 
2.80 
4.86 
5.50 
19.50 
3.93 
11.90 
2.45 
8.68 
5.15 
1.80 
3.27 
4.24 
1.29 
12.80 
2.54 
2.96 
8.57 
19.90 
20.00 
15.20 
8.64 
26.90 
8.52 
10.20 
20.50 
11.60 
9.98 
18.40 
24.90 
7.74 
19.20 
4.87 
17.00 
9.69 
3.38 
2.14 
38.70 
51.20 
28.90 
12.70 
42.70 
20.90 
20.40 
91.70 
2.23 
53.30 
24.10 
14.00 
9.92 
14.40 
19.00 
10.30 


KHCN/HNC) 


0.94 
1.34 
1.04 
1.20 
1.64 
1.16 
2.52 
1.54 
1.59 
0.52 
1.65 
0.50 
0.60 
1.99 
1.63 
0.47 
0.90 
0.97 
2.24 
1.21 
1.14 
2.46 
2.24 
1.23 
2.76 
1.01 
0.72 
2.35 
0.97 
1.83 
2.10 
2.24 
2.05 
1.55 
2.73 
2.10 
321 
4.29 
3.18 
1.17 
2.15 
1.70 
1.55 
0.58 
1.11 
1.66 
1.42 
2.01 
1.78 
0.56 
2.24 
2.14 
1.83 
2.28 
2.15 


а Integrated line intensity of the N3H* main component (1-0 Е = 2-1 F = 3-2). 
n Integrated line intensity of the C;H main component (1-0 3/2-1/2 F — 2-1). 


ІН?) 


8.74 
10.50 
11.20 
2.84 
6.11 
5.59 
13.90 
4.63 
13.20 
2.57 
6.51 
9.38 
5.94 
4.79 
4.78 
1.76 
2.11 
17.90 
19.70 
11.40 
14.90 
8.29 
25.50 
11.20 
13.20 
14.00 
12.50 
11.60 
22.50 
15.90 
3.34 
21.30 
2.51 
12.80 
6.91 
0.39 
19.70 
10.50 


22.50 
11.00 
35.00 


18.50 
7.67 
2.67 
TaT 


2.97 
4.06 


KHC3N) 


2.29 
3.49 
7.62 
0.22 
0.69 
0.97 
2.42 
0.74 
2.04 
0.40 
1.25 
1.76 
1.08 
0.34 
0.45 
1.86 
0.18 
9.75 
5.71 
4.43 
5.55 
1.13 
6.55 
2.41 
3.88 
5.03 
5.72 
3.22 
5.70 
5.68 
0.98 
4.56 
1.82 
3.99 
2.72 
0.35 
0.14 
15.30 
24.50 
11.20 
15.40 
27.20 
6.72 
6.78 
57.60 
19.00 
10.50 
4.35 
1.03 
3.13 
2.34 
1.64 


KHP?Cco^) 


2.01 
1.83 
2.46 
0.69 
2.01 
2.08 
2.63 
0.79 
2.60 
0.58 
0.85 
1.02 
0.95 
0.82 
0.40 
0.60 
0.47 
1.02 
4.80 
2.97 
3.50 
3.78 
1.24 
5.08 
2.20 
2.59 
4.85 
3.02 
3.20 
3.17 
3.65 
1.31 
2.55 
1.09 
3.26 
1.49 
0.43 
0.74 
8.91 
9.69 
6.10 
2.72 
14.50 
5.75 
4.51 
15.40 
6.49 
4.30 
3.59 
0.93 
2.46 
2.99 
1.28 


KCH)? 


1.52 
2.02 
3.14 
0.25 
1.55 
1.94 
5.12 
1.32 
2.11 
0.66 
1.34 
0.73 
0.89 
1.03 
1.07 
1.11 
2.51 
0.52 
1.20 
5.23 
2.87 
5.18 
4.28 
1.83 
7.14 
3.19 
3.93 
6.66 
2.95 
3.10 
6.11 
6.63 
2.56 
5.22 
3.07 
5.97 
4.13 
2.36 
1.28 
11.90 
12.60 
6.97 
6.14 
9.90 
5.06 
6.36 
21.8 
15.10 
7.42 
6.90 
2.61 
3.90 
5.42 
2.96 
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Table 4 

Kinetic Temperatures Obtained from the Different Methods (Units in K) 
Source Tx(HCN/HNC) Tease Tx (NH3) Source Тк(НСМ/НМС) Tasst Tx (NH3) 
IRDCOI 1.1 9.4 18.0 HMPO18566 9.7 22.8 19.9 
IRDCO28.1 13.4 17.0 HMPO19217 18.3 23.6 27.5 
IRDC028.2 10.4 17.0 HMPO19410 21.0 23.6 22.1 
IRDC048.6 12.0 19.0 EE НМРО20126 22.4 EE 23.1 
IRDCO079.1 16.4 15.0 15.3 HMPO20216 20.5 ee 19.2 
IRDCO079.3 11.6 14.0 HMPO20293 15.5 20.2 20.1 
IRDC18151 25.2 EE HMPO22134 27.3 24.5 20.9 
IRDC18182 15.3 20.0 HMPO23033 21.0 EE 18.5 
IRDC18223 15.9 18.0 HMPO23139 32.1 e. 27.7 
IRDC18306 3.2 19.0 HMPO23151 40.9 EE e. 
IRDC18308 16.5 19.0 HMPO23545 31.8 
IRDC18310 5.0 20.0 
IRDC18337 6.0 20.0 HMC009.62 11.7 25.1 
IRDC18385 19.9 20.0 НМСО10.47 21.5 23.2 te 
IRDC18437 16.3 20.0 HMC029.96 17 26.7 33.0 
IRDC184.1 47 20.0 HMC031.41 15.5 e. 20.3 
IRDC184.3 9.0 19.0 HMC034.26 5.8 28.4 er 
IRDC19175 9.7 20.0 HMC045.47 11.1 25.1 26.4 
IRDC20081 22.4 e HMC075.78 16.6 27:2 БЕ 
HMPO18089 12.1 23.4 23.0 UCH005.89 14.2 26.8 
HMPO18102 11.4 17.8 EE UCH010.10 20.1 21.8 
HMPO18151 24.6 ee 23.8 UCH010.30 17.8 26 
HMPO18182 22.4 21.5 te ЏСНО12.21 5.6 23.1 
HMPO18247 12.3 22.4 20.1 ЏСНО13.87 22.4 26.9 EE 
HMPO18264 27.5 20 22.2 UCH030.54 21.4 25.5 32.6 
HMPO18310 10.1 20 EE UCH035.20 18.3 EE 26.1 
HMPO18488 7.2 22 18.6 UCH045.12 22.8 31.3 БЕ 
HMPO18517 23.5 EE EE UCH045.45 21.5 29.3 25.4 


4.1. Testing the use of HCN/HNC Ratio as a 
Thermometer 


Employing simple molecules, ubiquitous in star-forming 
regions, to derive their physical parameters is a helpful strategy 
for methodically analyzing many regions and avoiding 
complex calculations. In this context, we use the HCN/HNC 
intensity ratio formulation proposed recently by Hacar et al. 
(2020) to estimate the kinetic temperature (Tk(HCN/HNC)). 
As done by the authors, we probed them through comparisons 
with dust temperature (Taus) and Тк obtained from the 
ammonia emission (7x(NH3)), when those parameters were 
available toward the sources of the analyzed sample. It is worth 
noting that the useful kinetic temperature range derived from 
the ammonia is between 15 and 40 K (Ho & Townes 1983; 
Urquhart et al. 2011), in coincidence with the valid range 
indicated for T«(HCN/HNC). 

Despite the limited amount of sources satisfying that the 
obtained Ty from the HCN/HNC ratio is >15 K with measured 
values of Тал: and T&(NH3) (see Table 4), some conclusion can 
be extracted. Also, from the sources with Тұ(НСМ/НМ№С) < 15 
K or lower than expected, we can obtain important conclusions 
about this tool. 


The comparison between Taust and Tk(NH3) (see upper panel 
of Figure 2) indicates that, within a range of 0-5 K of 
difference between both temperatures, the molecular gas and 
the dust could be considered that are thermally coupled, 
allowing us to suggest that in general this property could be 
extrapolated to the whole sample. This supports that the 
comparison between Tx(HCN/HNC), which is obtained from 
the gas emission, and the 74,4 can be made. 

Valid kinetic temperature values obtained following the 
HCN/HNC ratio (i.e., >15 K) seem to be more accurate for 
IRDCs and HMPOs than for HMCs and UCHII regions. By 
analyzing the average values (see Table 5), Ty obtained from 
the dust and ammonia emissions increases from IRDCs to 
UCHI regions as expected, but T,(HCN/HNC) increases only 
from IRDCs to HMPOs, and then it decreases in HMCs and 
UCHII regions. This suggests that in such sources, using the 
HCN-HNC tool to derive Ty can be not appropriate, and it may 
underestimate the actual temperature value. 

Hacar et al. (2020) calculated the Tx(HCN/HNC) toward a 
large sample of dense clumps extracted from the MALT90 survey 
(Foster et al. 2011; Jackson et al. 2013), and compared them with 
the Taust measurements derived by Guzmán et al. (2015). At this 
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Figure 2. Top-left: Kinetic temperature of ammonia derived from transitions (J, K ) = (2, 2)-(1, 1) versus dust temperature. Top-right: Dust temperature obtained from 
PPMAP output produced with the Hi-GAL maps (Marsh et al. 2017) versus kinetic temperature derived from the HCN-HNC integrated intensity ratio. Bottom: 
Kinetic temperature obtained from the ammonia versus kinetic temperature derived from the HCN-HNC integrated intensity ratio. The dashed black line, in all cases, 
indicates unity. Error bars in the T,(HCN/HNC) are displayed, and most of them are represented by the symbol size. 


point, they stated that the comparison between Тк and Taust should 
be done with caution given the fact that the beam resolution of the 
MALT90 targets and the measurement sensitivity within this 
beam do not allow resolving the temperature of individual clumps. 
Anyway, they pointed out that a correspondence between both 
temperatures prevails despite the mentioned observational caveats. 
However, they do not discuss the optical depth issue in the 
MALT90 sample as they did for their Orion results. This is an 
important matter to be considered, mainly when we work with a 
sample of different kinds of sources. As Hacar et al. (2020) 
mention, the increase of the HCN J — 1-0 line opacity, while the 
HNC J — 1-0 line likely remains optically thin, would reduce the 
HCN/HNC ratio, decreasing the values of Tk (HCN/HNC). This 
phenomenon could be occurring in the sources with 
Tx(HCN/HNC) < 15 K (see Section 4.1.1) and in some other 
sources that, although the Tę obtained from the HCN-HNC tool 
is greater than 15 K, their values are well below the temperatures 
obtained from the dust and/or ammonia (about 10 K of 
difference). This is the case of HMCs and UCHII regions in our 
sample, in which, in general, it is observed that the HCN 


Table 5 
Calculated Average Temperatures (Units in K) 

Source Taust Tx(NH3) Тк(НСМ/НМС) 
IRDC 18.5 + 0.4 15,34 18.5--Е 12 
НМРО 21.8 + 0.6 22.2 + 0.8 24.9 + 1.7 
HMC 26.0 + 0.6 26.6 + 3.7 16.3 + 1.6 
UCHII 26.3 + 1.1 28.0 + 2.3 19.8 + 1.0 
Note. 


а Value obtained from a single source. 


hyperfine line F — 1-1 appears absorbed, indicating high optical 
depths (for an example, see Figure A5). 


4.1.1. Sources with Ty(HCN/HNC) < 15 К 


In our sample, there are several sources such that the HCN/ 
HNC ratio yields kinetic temperatures lower than 15 K, which 
were neither included in the analysis nor the comparisons with 
Taust and Tk(NH3) performed above. These sources have 
I(HCN)/I(HNC) ratios close to unity or even lower than one, 
which according to Hacar et al. (2020), in those cases, the 
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Figure 3. Plots of the integrated line emission versus the 870 jum peak flux. Top panels display МН“ (main component, left) and СН (main component, right). 
Bottom panels feature: HNC (left) and HC3N (right). Error bars in the y-axis represent the formal error in the integration of the line emissions. 


uncertainties seem to increase, likely due to the combination of 
excitation and opacity effects. The percentages of sources with 
this issue are: 58%, 30%, 28% and 11% for IRDCs, HMPOs, 
HMCs and UCHII regions, respectively, suggesting that may be 
a condition that is more pronounced at earlier stages. 

Taking into consideration the HCN hyperfine line anomalies 
(Walmsley et al. 1982; Loughnane et al. 2012), which prevent 
obtaining reliable values of HCN opacities, we carefully 
analyzed each HCN and HNC spectrum to look for signatures 
of high optical depths. We found that the HCN spectra of 
sources with Tx(HCN/HNC) < 15 K have pronounced satur- 
ation and/or self-absorption features in some cases at both the 
main component and hyperfine lines (see Figure 7 for an 
example), and in other cases in the hyperfine line F = 1-1. In 
the case of the HNC, we did not find such spectral features. 

We point out that the use of the HCN/HNC ratio as a 
universal thermometer in the ISM should be taken with care. 
We suggest that such a thermometer could be used only in 
some IRDCs and HMPOs (when the derived Kinetic temper- 
ature is not lower than 15 K; see Section 4.1.1) and in more 
evolved regions, for instance, HMCs and UCHII regions, this 
tool underestimates the temperature. Thus, we conclude that the 


10 


HCN-HNC tool as a kinetic temperature estimator should be 
used only after a careful analysis of the HCN spectrum, 
checking that no line, either the main or the hyperfine ones, 
presents absorption features; otherwise, it is not an useful tool 
for calculating kinetic temperatures. 


4.2. An Exploration of Molecules as "Chemical Clocks" 


Molecular species such as HNC, HCN, HPCO*, С.Н and 
HC3N are relatively easy to be observed and they give us a 
significant amount of scope to search for differences in the 
chemistry as a function of the evolutionary stage in massive 
star-forming regions (Sanhueza et al. 2012; Yu & Wang 2015; 
Urquhart et al. 2019). After a similar analysis as presented by 
Yu & Wang (2015) performed on our sample of massive star- 
forming regions, in what follows we discuss our results. 

First, from the comparison between the molecular integrated 
intensities and the 870 um peak flux (see Figure 3), we found 
that, in general, the integrated intensities of N;H*, СН, HNC 
and HC3N increase with the submillimeter peak flux. In the 
case of IRDCs, HMPOs and HMCs, this increment correlates 
with the evolutionary stage of the sources, i.e., the position of 
each kind of source in the plots seems to be sectorized. UCHN 
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Table 6 
FWHM Av from Gaussian Fittings (Units in km s^!) 
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Source 


IRDCO11.1 
IRDC028.1 
IRDC028.2 
IRDC048.6 
IRDCO79. 1 
IRDCO079.3 
IRDC18151 
IRDC18182 
IRDC18223 
IRDC18306 
IRDC18308 
IRDC18310 
IRDC18337 
IRDC18385 
IRDC18437 
IRDC184.1 
IRDC184.3 
IRDC19175 
IRDC20081 
HMPO18089 
HMPO18102 
HMPO18151 
HMPO18182 
HMPO18247 
HMPO18264 
HMPO18310 
HMPO18488 
HMPO18517 
HMPO18566 
HMPO19217 
HMPO19410 
HMPO20126 
НМРО20216 
НМРО20293 
НМРО22134 
НМРО23033 
НМРО23139 
НМРО23151 
НМРО23545 
HMC009.62 
HMC010.47 
HMC029.96 
HMC031.41 
HMC034.26 
HMC045.47 
HMCO075.78 
CH005.89 
CH010.10 
CH010.30 
CH012.21 
CH013.87 
CH030.54 
CH035.20 
CH045.12 
CH045.45 


CIC Cc c e cu 


Notes. 


а Av of the N;H* main component (1-0 F, = 2-1 F = 3-2). 
> Av of the СН main component (1-0 3/2-1/2 F — 2-1). 


АУН?) 


3.03 
3.39 
2.25 
2.27 
2.30 
2.54 
3.08 
2.76 
3.41 
2.84 
3.82 
3.46 
3.10 
2.67 
2.88 
2.82 
1.91 
3.39 
5.00 
2.50 
3.55 
3.15 
3.31 
3.14 
3.48 
3.16 
3.94 
4.41 
2.50 
2.54 
2.20 
2.58 
2.17 
3.00 
2.83 
2.05 
4.41 


3.09 


4.41 
3.90 
3.50 
4.58 
3.14 
3.59 
3.29 


4.05 
3.20 


Av(HC3N) 


2.10 
2.87 
3.62 
1.31 
1.01 
1.20 
3.28 
2.17 
2.86 
2.43 
2.68 
2.97 
3.17 
1.74 
2.32 
2.52 


3.60 
3.67 
1.90 
3.05 
2.43 
3.14 
1.88 
3.20 
2.86 
3.62 
4.25 
1.85 
2.41 
1.33 
1.88 
1.25 
2.35 
2.50 
1.81 
1.34 
4.96 
7.48 
3:52 
5.99 
6.33 
3.93 
3.00 
3.94 
4.54 
7.87 
2.46 
2.42 
3.23 
3.98 
3.92 


Av(HPCO*) 


2.11 
2.69 
3.22 
1.52 
1.68 
1.68 
2.53 
1.61 
2.78 
1.74 
2.17 
2.11 
1.84 
1.30 
1.52 
2.53 
1.95 
1.16 
3.54 
3.65 
2.16 
3.02 
2.52 
2.93 
2.33 
3.15 
2.52 
3.73 
3.99 
1.80 
2.38 
1.37 
1.94 
1.69 
2.52 
2.48 
2.24 
2.34 
4.02 
7.06 
2.72 
4.90 
4.71 
4.31 
3.54 
3.64 
4.35 
5.94 
2.93 
2.55 
3.32 
3.39 
3.41 


Лу(СН)? 


2.11 
2.57 
3.53 
1.11 
1.79 
1.72 
3.85 
2.45 
2.87 
2.03 
3.00 
2.61 
2.85 
1.96 
1.88 
2.90 
4.60 
1.72 
1.56 
3.47 
3.49 
2.26 
3.29 
2.66 
3.01 
2.23 
4.61 
2.97 
4.20 
3.75 
1.99 
2.34 
1.82 
2.03 
1.84 
2.58 
2.92 
2.42 
1.80 
4.14 
6.95 
3.27 
7.03 
3.46 
4.33 
4.36 
3.33 
6.26 
6.14 
2.94 
3.18 
2.73 
4.00 
4.70 
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Table 7 
Average values of the Molecular Line Widths (Units in km 5 !) 
Source Av(N>H*) Av(HC3N) AwHPCO^) Av(C>H) 
IRDC 2.88 + 0.11 2.39 + 0.19 2.00 + 0.13 2.47 + 0.20 
HMPO 3.09 + 0.17 2.51 + 0.19 2.59 + 0.16 2.78 + 0.18 
HMC 3.95 + 0.31 5.02 + 0.62 4.46 + 0.51 4.79 + 0.58 


UCHI 3.62 + 0.20 4.04 + 0.60 3.64 + 0.38 4.19 + 0.49 


regions show the same trend, but their points are overlapped 
with those of the other sources in all cases except for the NH * 
case. Yu & Wang (2015) found that the NjH* integrated 
intensities of HII regions tend to be lower that those for 
MYSOs. This behavior seems to occur with the UCHII regions 
of our sample in comparison with the other sources, suggesting 
that the destruction path due to electronic recombination, where 
N5H* is destroyed by free electrons from the surroundings 
(Vigren et al. 2012), could be ongoing in the UCHII regions. 
This result complements what was presented by Yu & Wang 
(2015) with sources at earlier stages, and suggests that the 
relation between the molecular integrated intensities and the 
870 um peak flux could be useful to distinguish regions among 
IRDCs, HMPOs and HMCs. 

The comparisons between the column densities obtained 
from Gerner et al. (2014) also show an interesting trend 
regarding the evolutionary stage of the sources. Yu & Wang 
(2015) stated, based on astrochemical models (e.g., Bergin 
et al. 1997; Nomura & Millar 2004), that the НСО“ column 
density reflects the amount of H, density in a clump because its 
abundance seems to not vary much with time. In line with their 
work, we compared N(CH) and N(N5H^) with N(HP?CO") 
(see Figure 5). The N(HPCcO*) values of our sample range in a 
wider interval of values (from some 10? cm ? for IRDCs to 
values close to 10!^ ст 2 for HMCs) than the МН СО") 
presented by Yu & Wang (2015). An increase in the МСН) 
and N(N>H*) with the N(H'°CO*) is observed. In the case of 
N(C3H) versus МН СО") relation, such an increment also 
presents a conspicuous correlation with the evolutionary stage 
from IRDCs to HMCs in which the positions of each kind of 
source are sectorized in the plot. In the UCHII regions, both 
column densities have similar values compared with those of 
HMCs and the largest ones of HMPOs. A similar behavior, 
with larger dispersion, is observed in the comparison between 
ММН?) versus N(H'°CO*). The N;H* and H^CO* column 
density increment from IRDCs to HMCs could be explained by 
their progressive formation through the HZ mechanisms (H+ 
8СО — HPCO* + Hp, and Hj + № ^ N;H* + Н; Jørgen- 
sen et al. 2004). Then, the brake on the constant increment in 
the column density values in the UCHII regions compared with 
the previous sources can be due to the beginning of the 
destruction of such molecules by electronic recombination (Yu 
& Xu 2016). 
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Figure 4. Plots of the FWHM line widths (Av) obtained from Gaussian fittings to the emission lines of C,H, H'CO*, HC;N and N;H*. The dashed black lines 


indicate unity. Error bars represent the formal errors of the Gaussian fittings. 


It is worth noting that the discussed behavior in the column 
densities (which have molecular excitation assumptions) 
regarding the correlation with the evolutionary stage from 
IRDCs to HMCs is quite similar to the comparison between the 
integrated intensities and the peak submillimeter flux (which 
are direct measurements). Thus, in the line of a chemical clock 
analysis, even though more statistics are necessary, we propose 
that these relations, mainly the ММН?) versus NHPCO*) 
and the /(C;H) versus Езто ит. in which the different kinds of 
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sources (IRDCs, HMPOs and HMCs) are clearly separated, 
could be used to differentiate them, which can be useful in 
works handling a large amount of sources of unknown, or not 
completely known, nature. 

The comparison between the relative abundance ratios 
[C;H]/[H?CO'^] versus [NH *]/[H?CO''] (Figure 6) pre- 
sents a high dispersion, mainly along the y-axis, i.e., in the 
N(C>H)/N(H'°CO*) ratio, and does not show a remarkable 
difference as it was found between MYSOs and H II regions by 
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Figure 5. Column densities of CH (top) and NoH* (bottom) against column 
density of H^CO*. The column density values were obtained from Gerner 
et al. (2014). 


1e+14 
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Yu & Wang (2015). However, along the x-axis, i.e., in the 
N(N;H')/N(HPCO^) ratio, some trend can be appreciated, 
mainly such a ratio seems to be larger in IRDCs in comparison 
with the other kind of sources, which in the line of a chemical 
clock analysis, we suggest that this abundance ratio can be used 
to differentiate the earliest stage of the star-forming regions. 
The molecular line widths (FWHM) are usually associated 
with the kinematics of the molecular gas and can be affected by 
multiple events, such as shocks, outflows, rotation of the clump 
and turbulence. In this work, we assume that line widths 
originate mainly from the increasing turbulence that arises as a 
consequence of evolving star formation stages (Sanhueza et al. 
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Figure 6. Relative abundance of [C;H]/ (НСО versus [NH] / [H Cot]. 
The column density values were obtained from Gerner et al. (2014). 


HCN (1-0) — IRDC 011.1 


88.63 88.64 
Frequency [GHz] 


88.65 


88.62 


Figure 7. Example of an HCN spectrum with strong signatures of saturation 
which affects the use of the HCN-HNC ratio as a thermometer yielding 
unreliable temperatures below 15 K. 


2012; Yu & Wang 2015). With this in mind, by analyzing the 
average line width of МН”, HC4N, H'°CO* and СН (see 
Table 7), it can be appreciated that in the four species, the line 
width rises until reaching the HMC stage. Such an increase 
could be a consequence of the gas dynamics related to the star- 
forming processes that take place in the molecular clumps 
(Fontani et al. 2002; Pirogov et al. 2003). Particularly, our 
Av(N5H^) average values are in quite good agreement with 
those presented by Sanhueza et al. (2012): the Av(NH ^) 
average value that we obtained for IRDCs is similar to that 
obtained by the authors in their so-called "quiescent" and 
"intermediate" clumps, while the value that we obtained for 
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HMPOs is in agreement with the values of their “active” 
clumps. Finally, our Av(N>H*) average values obtained for 
HMCs and UCHII regions are slightly larger than the value 
obtained for their “red” clumps. Following Sanhueza et al. 
(2012), we can confirm that line widths of №Н” slightly 
increase with the evolution of the clumps. 

As found by Sanhueza et al. (2012) in their sample of clumps 
embedded in IRDCs, we observed that in our sources, the 
H'CO* and HC4N have slightly narrower line widths than 
N3H* (see panels d and f in Figure 4). This hints that they trace 
similar optically thin gas emanating from the internal layers of 
the regions, and according to our findings, it seems to be 
independent of the kind of source. Yu & Wang (2015) found 
that the best Av correlation is between C;H and N5H* for their 
sample of EGOs. In our case, we did not find such a good 
correlation between these molecular species (see panel e in 
Figure 4). The СН line widths are slightly narrower than those 
of the N;H *. Sanhueza et al. (2012) also remarked that the C,H 
appears to present the best correlation with NH * among all the 
studied Av relations. However, by inspecting their Av(CoH) 
versus Av(N>H*) plot, we observe a similar behavior as our 
plot, where most of the points tend to be below the unity line. 

In the spectra of these four molecular species toward the 
whole sample, we did not find line wings that may probe 
outflows, except for the HC4N spectra of HMPO18247, 
HMC029.96 (see Figure Аб) and ЏСН10.30, and in the 
H'CO* spectrum of HMC34.26, in which small spectral 
wings appear. UCH10.30 is an EGO; the other sources present 
4.5 um extended emission and diffuse emission at the Ks band. 
Besides these sources, it can be noticed that the near-infrared 
(near-IR) evidence of jets or/and outflows does not have any 
direct correlation with spectral features in the emission of 
МН”, HC3N, НСО? or СН. 


5. Summary and Concluding Remark 


We presented a spectroscopic molecular line analysis of 55 
high-mass star-forming regions with the aim of probing some 
chemical tools that were recently proposed to characterize such 
regions. The analyzed sources are classified as IRDCs, 
HMPOs, HMCs and ЏСНП regions, according to an evolu- 
tionary progression in the high-mass star formation. The main 
results can be summarized as follows: 

1. The emissions of the HCN and HNC isomers were used to 
estimate the kinetic temperature through a new “thermometer” 
whose formulation was proposed by Hacar et al. (2020), and we 
investigated its use in the presented sample of sources. By 
comparing the Ty derived from the HCN/HNC ratio with 
temperatures obtained from the dust and ammonia emission, we 
found that the use of such a ratio as a universal thermometer in 
the ISM should be taken with care. The HCN optical depth is a 
big issue to be taken into account. Line saturation may explain 
the derived temperature values below 15 K, which must be 
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discarded, and mainly correspond to IRDCs and HMPOs. In 
addition, we found that, although the 7, obtained from the 
HCN-HNC tool is greater than 15 K, its value could be far from 
the temperatures obtained from the dust or ammonia, yielding 
lower temperatures. This is the case of HMCs and UCHII regions 
of our sample, in which, in general, it is observed that the HCN 
hyperfine line F = 1-1 appears absorbed. In conclusion, we point 
out that the HCN-HNC tool as a kinetic temperature estimator 
should be used only after a careful analysis of the HCN 
spectrum, checking that no line, either the main or the hyperfine 
ones, presents absorption features. 

2. Additionally, we analyzed the molecular species NoH ^, 
HCN, НОСО“ and СН in the sample of sources focusing on 
the use of them as "chemical clocks." The comparison of the 
molecular integrated intensities of HC4N, МН”, C,H and also 
of HNC with the 870 um peak flux, as well as the relations of 
the column densities of C2H and МН" with that of the 
НСО", can be useful to distinguish regions among IRDCs, 
HMPOs and HMCs, which complement the Yu & Wang 
(2015) results. On the other hand, from the [C2H]/ [H's СОТ] 
versus [N9H *] / [HCO] relation, we can point out that even 
though the [C2H]/ [H?CO*] ratio shows a large dispersion 
along the analyzed sources, the [N;H*] / [H COH] ratio seems 
to be larger in IRDCs in comparison with the other kind of 
sources analyzed in this work. Finally, regarding the molecular 
line widths, we found that the Av rises from the IRDC to the 
HMC stage, which could be due to the increasing turbulence as 
a consequence of the evolution of the star-forming processes. 

In conclusion, this work explores chemical tools, some of 
them based on direct measurements, to be applied in different 
ISM environments. In that sense, we probed such tools in a 
new sample of sources with respect to previous works, and 
these results not only contribute to more statistics in the 
literature but also complement such works with other types of 
sources. Using direct tools, like, for instance, the ratios of 
different molecular parameters that can be directly measured 
from the observations, can be very useful mainly when a large 
sample of sources is handled. If it is proven that such tools are 
reliable, they can be used to obtain important statistical 
information in a simple way. Our work points to it and 
encourages performing similar works in larger samples of 
sources of different types. 


Acknowledgments 


We thank the anonymous referee for her/his useful 
comments that improved our work. N.C.M. is a doctoral 
fellow of CONICET, Argentina. S.P. is a member of the 
Carrera del Investigador Científico of CONICET, Argentina. 
This work was partially supported by the Argentina grants PIP 
2021 11220200100012 and РІСТ 2021-СКЕ-ТП-00061 
awarded by CONICET and ANPCYT. 


Research in Astronomy and Astrophysics, 24:015007 (23pp), 2024 January 


Appendix 
Example Spectra 


In this appendix we include some spectra of the analyzed 
molecular lines as an example of each kind of source, showing 
the integrated area in the case of the HCN and HNC emissions 
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and the Gaussian fittings for the H^CO^, СН, HC3N 
and N5H*. Spectra of IRDC18151 are shown in Figures Al 
and A2, of HMPO20126 in Figures A3 and A4, of 
НМСО.29.26 in Figures A5 and Аб, and of UCHI 45.45 in 
Figures A7 and A8. 
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Figure A1. IRAM spectra containing the HCN, НСО" and C;H lines toward IRDC18151. The integrated area in ће HCN emission and the Gaussian fittings for the 


other molecular lines are shown. 
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Figure A2. IRAM spectra containing the HNC, HC3N and N>H* lines toward IRDC18151. The integrated area in the HNC emission and the Gaussian fittings for the 
other molecular lines are shown. 
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Figure АЗ. IRAM spectra containing the HCN, НСО“ and CH lines toward HMPO20126. The integrated area in the HCN emission and the Gaussian fittings for 
the other molecular lines are shown. 
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Figure A4. IRAM spectra containing the HNC, HC3N and N2H* lines toward НМРО20126. The integrated area in the HNC emission and the Gaussian fittings for ће 
other molecular lines are shown. 
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Figure A5. IRAM spectra containing the HCN, H'*CO* and С.Н lines toward HMC029.96. The integrated area in the HCN emission and the Gaussian fittings for ће 
other molecular lines are shown. In the case of HCN it is remarked that the hyperfine component F = 1-1 is absorbed. 
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Figure Аб. IRAM spectra containing the HNC, HC3N and N5H* lines toward HMC029.96. The integrated area in the HNC emission and the Gaussian fittings for the 
other molecular lines are shown. In the case of the HC3N emission small spectral wings appear, and the Gaussian fitting corresponds to the main component. 
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Figure A7. IRAM spectra containing the HCN, НСО“ and СН lines toward UCHII 45.45. The integrated area in the HCN emission and the Gaussian fittings for 


the other molecular lines are shown. 
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Figure АЗ. IRAM spectra containing the HNC, HC3N and NoH* lines toward ЏСНП 45.45. The integrated area in the HNC emission and the Gaussian fittings for the 


other molecular lines are shown. 
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